Origin of non-exponential relaxation in a crystalline ionic conductor: a 
multi-dimensional ^^^Ag NMR study 
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The origin of the non-exponential relaxation of silver ions in the crystalline ion conductor 
AgyPaSii is analyzed by comparing appropriate two-time and three-time ^"^Ag NMR correlation 
functions. The non-exponentiality is due to a rate distribution, i.e., dynamic heterogeneities, rather 
than to an intrinsic non-exponentiality. Thus, the data give no evidence for the relevance of corre- 
lated back-and-forth jumps on the timescale of the silver relaxation. 
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Solid ion conductors are interesting materials from 
the viewpoints of both fundamental and applied science. 
While fast ion transport has important applications in 
energy and information technologies, many of its dy- 
namic aspects are poorly understood at the atomic level. 
A prominent feature of many solid ion conductors is their 
non-exponential relaxation behavior, reflecting complex 
ion dynamics P, 0, ■ Two fundamentally different ex- 
planations are possible 0. In the heterogeneous sce- 
nario, all particles are random- walkers, but a distribu- 
tion of correlation times G(lgr) exists. In the homo- 
geneous scenario, all particles obey the same relaxation 
function, which is, however, intrinsically non-exponential 
due to correlated back-and-forth jumps. For solid ion 
conductors, the concept of correlated forward-backward 
motion is one of the cornerstones of current modelling 
approaches. It is suggested by a strong frequeiicy depen- 
dence of the electric conductivity a{u) H, S 0- Thus, 
it might be tempting to conclude that the homogeneous 
scenario applies to ion dynamics, i.e., ion transport might 
be slowed down by extended back-and-forth jump se- 
quences. 

Despite the technological importance of solid ion con- 
ductors, the origin of the non-exponential relaxation has 
remained unresolved to the present date. This is due 
to the fact that two-time correlation functions (2T-CF) 
such as scattering functions or conductivities are intrinsi- 
cally unable to decide between both scenarios, but rather 
three-time correlation functions (3T-CF) are required 0. 
The lifetimes of dynamic heterogeneities, if they exist, 
can be measured by four-time correlation functions (4T- 
CF). As was shown for applications in supercooled liq- 
uids 0,0,0, all of these multi-time correlation functions 
(MT-CF) are available from multi-dimensional NMR ex- 
periments. There, the 3T-CF has a simple interpretation: 
a dynamic filter meant to select slow particles is applied 
in the first time interval and the relaxation of the selected 
subset is observed in the second time interval. If dynamic 
heterogeneities exist, such a selection is possible and the 



3T-CF decays more slowly than the regular 2T-CF. 

For solid ion conductors, previous 2T-CF NMR stud- 
ies have confirmed that relaxation is non-exponential 
[ifil Ull IT3. [l3j |. Moreover, it has been shown that 
MT-CF monitoring the diffusion of single silver ions can 
be recorded in ^''^Ag NMR In this letter, we use 

multi-dimensional ^'^^Ag NMR to study for the first time 
the nature of the non-exponential relaxation in a solid 
ion conductor in detail. Specifically, we measure 3T- 
CF in polycrystalline AgyPaSn to resolve the homoge- 
neous and the heterogeneous contributions to the non- 
exponential silver ionic relaxation. Moreover, the lifetime 
of the dynamic heterogeneities is quantified by means 
of 4T-CF. 7-Ag7P3Sii exhibits a high dc conductivity 
adc = 2.5 X lO-^'S/cm at T « 300K jlj. We study sil- 
ver dynamics in the /3-phase below a phase transition 
at T = 209K To improve the signal-to- noise ratio 

Ag7P3Sii was synthesized starting from ^°^Ag enriched 
(97%) metal powder. Further experimental details are 
described in Refs. [ilE3- 

In soUd-state ^^^Ag NMR, the chemical shift (CS) in- 
teraction dominates the observed frequency shift, u. The 
CS tensor describes the magnetic shielding of the applied 
static magnetic field Bq at the nuclear site due to neigh- 
boring electrons. In our case, fast silver rattling motions 
in the potential minimum affect the CS tensor so that a 
single, time-independent to can be ascribed to each sil- 
ver site. While the different local environments result in 
distinguishable values of w, the NMR frequencies at sites 
related by translational symmetry are identical. Thus, 
any change in the resonance frequency is due to silver 
jumps between distinguishable sites. In other words, the 
sites within the unit cell provide a small number of NMR 
frequencies which are discontinuously adopted by the sil- 
ver ions. i°^Ag NMR MT-CF directly probe the time 
dependence of uj and, hence, silver diffusion on the ms-s 
timescale. Since the values of oj associated with the sil- 
ver sites depend on the orientation of the crystal with 
respect to Bq a powder average is observed for our poly- 
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crystalline sample. 

The measurement of NMR MT-CF is well described in 
the literature [1, 0] . In general, suitable multi-pulse se- 
quences are applied to manipulate the spin system. The 
pulses of these sequences divide the experimental time 
into a series of alternating short evolution times ip <Ct 
during which the respective NMR frequencies are de- 
tected and long mixing times tm ~ t during which dy- 
namics may take place. The lengths and the phases of 
the pulses depend on the spin / of the observed nucleus. 
In our case / = 1/2, the multi-pulse sequences can be 
adopted from NMR experiments [111113. In detail, 
the three-pulse sequence Pi — tp — P2 — tm — P3 — tp, or 
stimulated-echo sequence, is used to record ^^^Ag NMR 
2T-CF. Varying tm for constant tp we measure 

Fi'^itm) ( cos[cjitp] cos[w2^p] ) and 
F2'^{tm) oc (sin[Li;iip] sin[a;2ip] )• 

Here, loi and L02 are the frequencies during the two evolu- 
tion times separated by the mixing time tm ■ The brackets 
(. . .) denote ensemble averages. Suitable seven-pulse se- 
quences can be applied to correlate the frequencies wi, 
UJ2, UJ3 and UJ4 during four evolution times tr, separated 
by three mixing times t„ii, im2 and tms 
record the 4T-CF 
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F4,{tml,trn2, t-ms) OC ( COs[(tJ2-tJl)tp ] C0s(tJ3tp) COii{uJitp) ), 

which has proven very useful for a study of the lifetime 
of dynamic heterogeneities. A 3T-CF suited to analyze 
the origin of the non-exponential relaxation is obtained 
for tm2^0 in the seven-pulse sequences, i.e., ijj3=ijj2: 

F3{trnl,t„i3) OC ( COS [(cJ2 - t^l )ip ] C0s{uJ2tp) C0s{uj4tp) ). 

A straightforward interpretation of these MT-CF is 
possible when each silver jump to a distinguishable site 
leads to a complete loss of correlation. This condition 
can be met when large evolution times are applied so that 
Acutp ^ 2tt is fulfilled for basically all frequency changes 
resulting from these jumps. For Acutp ^ 2tt, it is easily 
seen using the trigonometric addition theorems that 

F2{tra) = F^'{tm) = F^'{tm) K ( C0s[(lJ2 "^l)ip] ) (1) 

where the term cos[{Lj2~iOi)tp] approximates the 5 func- 
tion 5{uj2-i^i) 8, 9]. In i"9Ag NMR studies of AgyPaSn, 
Eq. mis fulfilled for tp > lOO/zs [ij. For such evolution 
times, independent of the actual value of tp, -Fjl^m) quan- 
tifies the fraction of ions that occupy the same or - due 
to the identical u - any periodic site after a time tm- 
Similarly, one can rewrite 

F3{tmhtm3)c^{cos[{uj2-u}i)tp]cos[{uj4,-uj2)tp]) and (2) 

F4{tmhtm2,tm3) {C0s[{uj2-U}l)tp] COs[{uj4-UJ3)tp]) . (3) 

In studies of glassy ion conductors, i.e., in the absence of 
periodic sites, the term cos[(ti;2 — wi)tp] acts as a perfect 
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FIG. 1: F2"{tm) for polycrystalline AgyPaSn at various tem- 
peratures (tp = 100/j,s). We compare results for samples con- 
taining 48% ^°^Ag (solid symbols, Ref. [ll|) and 97% ^°^Ag 
(open symbols), respectively. The dashed lines are interpola- 
tions with a modified KWW function. Inset: Mean correla- 
tion times < T > together with an Arrhenius interpolation. 



dynamic filter which selects immobile ions during tmi 
{bj2 — ljJi)- Then, the relaxation of these selected ions is 
observed via F3{tm3) so that F3{tm3) ^ F2{tm3) directly 
indicates the existence of dynamic heterogeneities. Fur- 
ther, the lifetime of dynamic heterogeneities can be mea- 
sured when Fi{tm2) for tmi — ims = const, is recorded. 
In this experiment, identical filters are applied to probe 
the dynamic states of an ion during tmi and a time tm2 
later during tm3 so that an exchange of the dynamic state 
results in a decay of Fi{tm2), cf. below. Here, we will 
demonstrate that although the lattice periodicity leads 
to imperfect dynamic filters in applications on crystalline 
ion conductors i^3(tm3) and -F4(im2) still yield valuable 
insights into the nature of silver dynamics in AgyPaSn. 

In Fig. 1, we show ^o^Ag NMR 2T-CF for AgyPaSn 
at selected temperatures. The silver jumps lead to tem- 
perature dependent, non-exponential decays of F2{tm)- 
Experiments for the labelled sample and for a sample 
with natural isotopic abundance yield identical re- 
sults. Fitting to a modified Kohlrausch- Williams- Watts 
(KWW) function [l^, (1 - C) exp[-(t/r)^]-H C, the non- 
cxponcntiality is reflected by a stretching parameter of 
/3 = 0.42. The plateau value C « 0.12 « 1/8 indicates 
eight occupied, magnetically distinguishable silver sites 
UJ. The mean correlation time can be calculated ac- 
cording to < r > = (t//3) r(l//3) where r(a::) is the F 
function. Its temperature dependence is well described 
by an Arrhenius law with activation energy Ea = 0.37eV, 
see inset. We carefully checked that neither spin-lattice 
relaxation (Ti w 30s) nor spin diffusion affect the data for 
mixing times shorter than about a few seconds. 

The origin of the non-exponential silver dynamics in 
Ag7P3Sii is revealed by F3{tm3)- It is shown in Fig. 2 for 
the dynamic filters tmi = 10ms, 40ms and T — 176K. In 
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A t„=10ms 



-F3 {tmh 



:l/2(/(<„a+im3)+/(tml)/(im3)). (5) 
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FIG. 2: i^3(tm3;tmi) for AgrPaSii at T = 176K (tp = 
200/is). The results for tmi = 10ms, 40ms are compared with 
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0) ~ F2''{t,„-i) ^ F2{tm3)- The respective ex- 



pectations for purely homogeneous dynamics (dashed lines) 
and purely heterogeneous dynamics (solid lines) are included. 
The expectations for heterogeneous dynamics in the O-model 
(black) and the D-model (gray) are nearly identical. 



Here, the brackets denote the average over the rate dis- 
tribution G(fc) that is present in the case of dynamic 
heterogeneities, but becomes a delta-function for purely 
homogeneous dynamics. 

The homogeneous limit can be discussed in a model- 
free way. Since a rate distribution is absent, one directly 
has p'^ (i) = F2 (t) and, thus, 

F^"^"" (tmh tjns) — (l/2)[-F2(iml+^m3) + -p2 (t„il )^2 (im3)] ■ 

(6) 

For the heterogeneous limit, we first discuss the O-model. 
In this case, one can specify p'^(i) from simple argu- 
ments. For a jump rate k and one-dimensional dynamics, 
the probability to be at a site with the same NMR fre- 
quency is (1/2)[1 -l-expf— 2fct)]|l^. Since the dynamics 
along the different dimensions are independent of each 
other p'=(i) = {(l/2)[l-|-exp(-2fci)]}3 follows. Using this 
expression for p'^{t) one can determine the distribution 
G(fc) from Eq. ^ and use it to calculate p^'^^ '-' via Eq. 
El Now, we turn to the heterogeneous scenario of the 
D-model. For a single jump rate fc, one has 
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exp(-fci„). 



(7) 



particular for the longer tmi: F-i{tmd) decays much more 
slowly than i^2(im3)- In studies of supercooled liquids, 
the theoretical expectations for purely heterogeneous dy- 
namics, i^3'^*(im3), and purely homogeneous dynamics, 
F^°™(tm3), can be expressed in a model-free way based 
on the measured 2T-CF [llESli^- For Ag7P3Sii, the 
presence of a small number of discrete NMR frequencies, 
which are discontinuously adopted, complicates the anal- 
ysis. In general, one has to specify the crossover paths 
between the different sites/cj to define the purely hetero- 
geneous and the purely homogeneous scenario. Since de- 
tailed information is not available, we restrict ourselves to 
the analysis of two extreme cases. For reasons of simplic- 
ity, we assume that there are eight distinguishable sites 
on a cubic lattice. While the frequencies of every other 
site in each dimension are identical in the ordered (0-) 
model, the iV w 8 values of uj are randomly distributed in 
the disordered (D-) model. 

For the following calculations, it is useful to rewrite 
F3 a l/2{(cos[(cji-Li;4)tp] ) -I- {cos[{ijJi+ ijJA-2ui2)tp] )}. 
Since Aa;ip3>27r, ions showing lui~lu4 and cji =0^2 = t^4 
contribute to the first and the second term, respectively. 
We denote p'^it) as the probability that after a time t 
an ion occupies a site with identical NMR frequency. If 
dynamic heterogeneities exist, this probability depends 
on the jump rate as indicated by the index k. In general, 
we can write 



With this choice oip''{tm), Eq. [Sjcan be rewritten as 
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[tml, tmz) — F2(tmi+ tma) + W [-^2 (irjxl ) + 



F2{tm) = {P^{tm)) and 



(4) 



F2{tm3) — F2{tml+ Unz) — '^]- (8) 

Actually, we use A^ = 8.4 for the D-model where this value 
results from a fit of Eq. ^ to the experimental data. 

In Fig. 2, we see that the experimental data for t,ni = 
10ms, 40ms agree well with F^'^*(im3) obtained within 
both the O-model and the D-model. The minor devia- 
tions in the plateau regime (im3 > Is) are likely due to an 
onset of spin diffusion. In contrast, F!}°"^{tm2,) distinctly 
deviates from the experimental results. The agreement 
of both limiting cases with respect to the heterogeneous 
limit clearly shows that, as is the case for supercooled 
liquids, the heterogeneous limit can be characterized in 
a basically model-free manner. We conclude that corre- 
lated back-and-forth jumps on the timescale of the silver 
relaxation are not a relevant feature of the dynamics. 

Finally, the lifetime of the dynamic heterogeneities in 
AgrPsSii is measured by recording F4^{tm2) for tmi = 
tm3 = ts ~ <T >. If there were perfect dynamic filters, 
slow silver ions (t > tg) would be selected during tmi 
and a time t„i2 later it would be checked whether these 
ions are in the same dynamical state during t^s, cf. Eq.O 
Thus, F4(tm2) would decrease when initially slow ions be- 
come fast (r <is) during tm2 until the re-equilibration of 
the respective subensemble is complete. In our case, the 
signal results from all ions showing cui = L02 and = a'4 
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FIG. 3: F4{ts,tm2,ts) for AgrPsSii at r = 176K {tp^200fj.s). 
The filter times ts are indicated. Lines: Interpolations with 
a modified KWW function. The inset shows the scaled mean 
decay times < rp^ > / < tfs >■ For comparison, one data point 
for r=182K is included (X). For ts = 2ms and for r = 182K 
< TF3 > has been estimated from Fg'^^, as calculated using 
Eq.lHl 



so that fast ions that occupy periodic sites at the rele- 
vant times contribute, too. Detailed random-walk sim- 
ulations based on the O- and the D-model have shown, 
however, that despite the imperfection of the dynamic 
filter Fi{tm2) is well suited to measure the timescale of 
exchange processes between fast and slow ions. 

The curves F4(tm2) for AgrPaSn at r= 176K and var- 
ious tg are shown in Fig. 3. Non-exponential decays to 
different plateau values for tm2 ~^ oo indicate exchange 
processes between slow and fast silver ions until the re- 
equilibration of the respective subensemble is complete. 
To quantify the exchange rate the decays are fitted to 
a modified KWW function. We obtain mean time con- 
stants <TFi > =390— 660ms that do not significantly de- 
pend on the filter time ts- These values are on the same 
order as < Tp^ > and, hence, the exchange processes be- 
tween slow and fast ions of the distribution occur on the 
timescale of the slow silver jumps. However, as indicated 
by < tf4 > / < tf3 > « 6, see inset, it takes a few rather 
than a single jump relaxation process to loose memory 
about the mobility, suggesting that high and low energy 
barriers are not randomly distributed in the unit cell of 
AgrPsSii. 

In summary, the observation of NMR MT-CF pro- 
vides valuable insights into the mechanism of ion mo- 
tion, not available by other experimental methods. The 
silver jumps in the crystalline ion conductor AgyPaSn 
are governed by a broad rate distribution within which 
exchange processes between slow and fast ions take place 
on the timescale of the slow jumps. Further, the purely 
heterogeneous scenario applies to the relaxation of the 



non-fast silver ions and, hence, back-and-forth jumps are 
not relevant on the time scale of the depopulation of the 
silver sites. At first sight, this seems to contradict the 
presence of back-and-forth motions inferred from the fre- 
quency dependent electric conductivity a^v). Both re- 
sults might be reconciled if one assumes that only very 
fast ions contribute to the dispersion of C7{v), whereas 
ionic motion on the typical timescale of the silver relax- 
ation is random- walk like. This conclusion is backed up 
by MD simulations on glassy ion conductors where the 
back-jump probability was found to depend strongly on 
the waiting times at the sites 0, . 
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through the Emmy Noether-Programm. 



[9 

[lo; 

[11 
[12; 

[is: 

[14 

[is: 
[16: 
[17: 

[is: 

[19 
[20 



C. T. Moynihan, L. P. Boesch and N. L. Laberge, Phys. 
Chem. Glasses 14, 122 (1973) 

C. Liu and C. A. AngeU, J. Non-Cryst. Sohds 83, 162 
(1986) 

D. L. Sidebottom, P. F. Green and R. K. Brow, J. Non- 
Cryst. Solids 183, 151 (1995) 

R. Bohmer, R. V. Chamberlin, G. Diezemann, B. Geil, A. 
Heuer, G. Hinze, S. C. Kuebler, R. Richert, B. Schiener, 
H. Sillescu, H. W. Spiess, U. Tracht and M. Wilhelm, J. 
Non-Cryst. Solids 235-237, 1 (1998) 

A. K. Jonscher, Nature 267, 673 (1977) 

M. D. Ingram, Phys. Chem. Glasses 28, 215 (1987) 

B. Roling, A. Happe, K. Funke and M. D. Ingram, Phys. 
Rev. Lett. 78, 2160 (1997) 

K. Schmidt-Rohr and H. W. Spiess, Multidimensional 
Solid-State NMR and Polymers, Academic Press, London 
(1994) 

R. Bohmer, G. Hinze, G. Diezemann and E. Rossler, 
Prog. Nucl. Magn. Reson. Spectrosc. 39, 191 (2001) 
R. Bohmer, T. Jorg, F. Qi and A. Titze, Chem. Phys. 
Lett. 316, 419 (2000) 

M. Vogel, C. Brinkmann, H. Eckert and A. Heuer, J. 
Non-Cryst. Solids 307-310, 971 (2002) 
M. Vogel, C. Brinkmann, H. Eckert and A. Heuer, Phys. 
Chem. Chem. Phys. 4, 3237 (2002) 

F. Qi, T. Jorg and R. Bohmer, Solid State Nucl. Magn. 
Reson. 22, 484 (2002) 

Z. Zhang and J. H. Kennedy, J. Elektrochem. Soc. 140, 
2384 (1993) 

C. Brinkmann, M. Vogel, H. Eckert, A. Heuer and A. 
Pfitzner, unpublished results 

U. Tracht, M. Wilhelm, A. Heuer and H. W. Spiess, J. 
Magn. Reson. 140, 460 (1999) 

R. Kohlrausch, Pogg. Ann. Phys. 4, 56 (1854); G. 
Wilhams and B.C. Watts, Trans. Faraday Soc. 66, 80 
(1970) 

A. Heuer, U. Tracht, S. C. Kuebler and H. W. Spiess, J. 
Molec. Structure 479, 251 (1999) 

U. Tracht, PhD Thesis, University of Mainz, Germany, 
1998 

A. Heuer and K. Okun, J. Chem. Phys. 106, 6176 (1997) 



5 



[21] H. Lammcrt, M. Kunow and A. Heuer, Phys. Rev. Lett. [22] M. Vogel, Phys. Rev. B (in press) 

90, 215901 (2003) 



